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Phases of the two-band model of spinless fermions in one dimension 
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We study the two-band model of spinless fermions in one dimension for weak repulsive interac- 
tions. In this case, the model is equivalent to the weakly interacting spinless two-leg ladder. We 
obtain non-universal analytic expressions for the power-law decays ex x ' of the charge-density- and 
oc a; -1 ' 7 of the superconducting pairing correlation functions. Leading order in the doping away 
from half filling 8 and tj_/t we find 7 = 1 + (it 2 /8)(t±/t) 2 5 2 (t and t± are the hopping terms along- 
and between the chains), such that superconducting pairing correlations dominate. We furthermore 
show that the transition from the superconducting phase to the usual one-dimensional (Luttinger) 
metal occurs via a mixed phase, where superconducting pairs are formed in the bonding band only. 
We give the phase diagram as a function of temperature and doping. 

PACS numbers: 71.10.Pm, 74.20.Mn 
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I. INTRODUCTION 

One-dimensional (ID) electron systems have attracted 
much attention over the past decade (s). While 3D sys- 
tems are well described by the Fermi liquid theory, for 
ID systems this is not the case. The low-energy physics 
is rather the one of a Luttinger liquid (LL) GL0], where 
the excitations are collective zero sound modes. The 
Hubbard- and the t — J chains provide so far the most 
reasonable description of ID (super)conductors [§Q. 

The extension to quasi-one-dimensional systems is ob- 
tained by coupling N (Hubbard, t — J) chains to form 
a AT-leg ladder. The interest in doped ladders started, 
when the possibility of superconductivity in two-leg lad- 
der materials was proposed [g-[7]] ■ For the Hubbard- and 
the t—J two-leg ladders different numerical methods have 
shown that ci-wave like superconducting correlations in- 
deed dominate |q-fiO||. Experimentally, superconductiv- 
ity has up to now only been observed in a two-leg ladder 
material under high pressure |llf| . 

The Hubbard chain is exactly solvable for all values 
of the hopping parameter t and the on-site repulsion U 
|p.2| , and the t — J chain at the supersymmetric point, 
t = J p3[ . Analytic works for ladders were made for 
particular limits: For weakly coupled chains (small in- 
terchain hopping t± <C U,t) 1 14 -16 1, and for stronger 
coupled chains, but small interactions (small on-site re- 
pulsion U <C t±,t, such that the AT-leg ladder is equiv- 
alent to a TV-band model) |14 17|-pi|. Finally, at- and 
very close to half filling, where recent works have shown 
how to map the (weakly interacting) two-leg ladder on 
the Gross-Neveu model [|19[ (for an application, see Ref. 
p(J) and the t—J—J± model on the XXZ spin chain pl| . 

In order to determine whether the system exhibits su- 
perconductivity or not, it is necessary to calculate the 
fluctuations (correlation functions) for the charge-, and 
spin-density and for the superconducting (SC) pairing. 
In the two-leg ladder, the decay of the charge density 



correlation function is oc x~ 2K and of the pairing cor- 
relation function oc x^ 1 ' ( 2K >, where K is the Luttinger 
liquid parameter (LLP) being one in the noninteracting 
case, K = 1. For sufficiently weak interactions pairing 
correlations thus dominate which is interpreted as the 
appearance of superconductivity p6| , p^| . 

Spinless fermions are usually studied as a model for 
the physical more interesting (but more difficult) case of 
spin-1/2 fermions — in particular when focusing on the 
interchain hopping t± 22 or the metal- insulator tran- 
sition |p3[ . Physically, spinless fermions can be consid- 
ered as completely polarized spin-1/2 fermions in a (high) 
magnetic field. For weak interactions, a two-leg ladder of 
spinless fermions is conveniently mapped on a two-band 
model. Without making the link to ladders, such a model 
has been treated in Ref. p4| . More recent works have 
studied two weakly coupled chains, ti < f |1J,G3,E5|. 
The authors found that the LL fixed point of the two 
separated chains is unstable upon weak coupling. But in 
contrast to the spin-1/2 case, in a two-leg ladder of spin- 
less fermions, the decay of the charge density correlation 
function is oc x~ K and that of the SC pairing is oc x~ x l K , 
and superconductivity does not occur for purely repulsive 
interactions (in the t — V—U model V > and U > 0) 
andi ± << pira. 



In this paper, we study the low-energy physics of the 
spinless two-leg ladder for weak repulsive interactions and 
finite interchain hopping t±_ ~ t. We obtain a non- 
universal analytic exponent for the power-law decays 
oc a;~ 7 of the charge density- and oc a; -1 / 7 of the SC 
pairing correlation functions. The phases as function of 
the temperature T and the doping away from half fill- 
ing S (for spinless fermions, < 5 < 0.5) are shown in 
Fig. 1. At half filling the ladder is insulating; upon dop- 
ing away from half filling, we come to a phase with one 
gapless mode, which exhibits due to the finite t± domi- 
nant SC correlations with mierchain pairing. When the 
doping is increased, this phase undergoes a transition to 



a phase with two gapless modes, where first SC correla- 
tions coexist with charge density correlations. The SC 
pairing correlations finally disappear completely at high 
doping, where the ladder becomes a ID Luttinger metal. 
We thus have a nontrivial crossover from the SC state to 
the normal conducting one. 
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FIG. 1. Phase diagram of the spinless two-leg ladder for 
t±/t = I and U/t = 0.2. The doping away from half filling 
is 5 and T is the temperature (the To is a high temperature 
cutoff). The solid line shows the crossover to a supercon- 
ducting (SC) phase with inierchain pairing, and the dashed 
line displays where the crossover from a metallic phase with 
charge density wave (CDW) excitations to a mixed phase with 
coexistence of SC and CDW excitations occurs. The dotted 
line separates the region where both bands are partially filled, 
from the region where one band is empty. At half filling, the 
ladder is insulating. 

To determine the low-energy properties, we make use 
of the renormalization group (RG) method [g7j and of 
bosonization pp^j. 

In section II, we introduce the Hamiltonian and prove 
the existence of a gapless mode using a modified version 
of the Lieb-Schultz-Mattis (LSM) theorem j^-JH). In 
section III, we give the RG-equations (RGEs) and the 
scaling of the couplings. In section IV, using bosoniza- 
tion the phases are derived and discussed. Appendix A 
contains the analytic solution of the one- loop RGEs. In 
Appendix B bosonization is explained and in Appendix C 
the correlation functions are calculated. 



II. THE TWO-LEG LADDER 

The noninteracting two-leg ladder of spinless fermions 
is given by the Bamiltonian, 

Ho = -t^2d](x + l)dj(x) + B.c. 

-t ± ^4(x)-i 2 (x)+B.c, (1) 



where t and t± are the hopping amplitudes along- and 
between the chains and dAx) creates a fermion in the 
chain j at the rung x. We are going to consider small 
repulsive interactions < U <C t, £j_. In this limit, it is 
a good approach to diagonalize first Hq by a canonical 
transformation, 



1 



^j(x) = -^(dx(x)±d 2 (x)). 



(2) 



Going over to the momentum space, we find a decoupling 
into two bands (we set the lattice parameter equal one), 



where the dispersion relations are 

tj(k) = =f£j_ - 2tcos(k). 



(3) 



(4) 



Band 1 is the bonding- and band 2 the antibonding 
band. By analogy with the 2D case, the associated trans- 
verse momenta are denoted as k± = 0, n. Since we are 
discussing only the low-energy physics, we linearize Cj 
around the Fermi momenta ifep-j resulting in Fermi ve- 
locities Vj — 2ts'm(kFj). For the operator *$>j at ±kFj, 
we write ^R/Lj{k) = ^j(±kFj + k). 

In the small U limit, it is a good approximation 
to take for kpj the values obtained in the noninter- 
acting system |17J]. The definitions of the band fill- 
ing n = (kpi + fcp 2 )/(27r) and the chemical potential, 
ji = ei(fcpi) = £2(^2) allow to calculate kpj and Vj as a 
function of n, t, and t± , 



i>i 2 = 2tsin 



nn ± arcsin 



t± 



2tsiii(jrn) 



(5) 



For spinless ladders, the band filling is < n < 1 and 
the hole doping away from half filling is S = 0.5 — n. The 
effect of the interchain hopping t± is thus included in the 
velocities Vj [0 . We like to point out that the velocities 
are not equal for finite ij_, 



2t 



Vi — V2 



tan(7m) 



(G) 



This is different to previous treatments of the spin- 
less two- leg ladder 22 . We will see that this differ- 



ence in the velocities has the remarkable effect of driv- 
ing the system to a SC state for repulsive interactions, 
U > 0. We note that both bands are partially filled, 
when t± < 2tsm(irn) 2 . 

We do not consider the half-filled case, kp X + kp 2 = it 
(vi = V2), where the ladder is insulating, allowing us 
to neglect umklapp processes (we also exclude the par- 
ticular points kpj = 7r/2). Including all interactions 
allowed by symmetry (leaving away completely chiral 



one's), in momentum space, the Hamiltonian is given by 
H = H + Hint, where 



Ho 



and 



i=i,2 



dkk 



*i,-(*)*iw(*) - *L(*)*«(*) 



(7) 



fca - fc? 



Hint = / dk\dk 2 dkzdk4 i 8{k 

ffl*Jjl( fc l)*-Rl( fe 2)*i 1 (fe3)*il(fc4) +.92(1 «-> 2) 
+3* (*L( fc l)^l( fc 2)*L( fc 3)*L2(fc 4 ) + 1^2 
+3* (*L( fc l)^2(fc 2 )*L( fc 3)*L2(fc 4 ) + 1^2 



(8) 



The bare values of the couplings are chosen as g\ = g 2 = 
f]x = 9t = U > 0. We will see, that the <?t interaction 
(pair hopping of left/right going quasiparticles from band 
one to band two) is the most relevant in determining the 
low-energy physics. 

Next, we apply the generalized LSM theorem (see Refs. 
p9|-pl[), to the Hamiltonian H in order to show the ex- 
istence of gapless modes. The particular symmetry of 
the Hamiltonian H allows to define "twist operators" for 
left/right going fermions separately, i.e., 



U R = exp 2m^2j^ Rj (x)^ Rj (x) 



x,3 



(9) 



where L is the length of the ladder (and similarly Ul 
for the left going fermions &lj)- Commuting U R / L 
with the translation operator T in real space, we obtain 
TUrT' 1 = e- 4 ™ R U R . The chiral symmetry n = 2v R = 
2vl proves then the existence of gapless (charge density) 
excitations at a wavevector 2irn = kpi + kp 2 , if n is not 
an integer (remember that we left away umklapp interac- 
tions in H present at half filling; they break the individ- 
ual Ur/l "quasisymmetry" ) pl[ . We show, that this is 
in agreement with the results obtained by bosonization 
(see section IV). However, the LSM theorem does neither 
tell us the number of gapless modes nor whether finally 
charge density- or superconducting excitations dominate. 



III. THE RG-EQUATIONS 

We give the RGEs resulting from the Hamiltonian H 
and the flow of the couplings depending on the ratio of 
the velocities V\/v 2 . We find that for Vi/v 2 < 7, all cou- 
plings diverge, while for V\/v 2 > 7, g t — > and the rest 
remains of the order of U. 

The model with the couplings g\, g 2 , and g x alone is 
exactly solvable (by bosonization, see below) and in par- 
ticular at a RG fixed point. Products of couplings without 



at least one gt do therefore not appear in the RGEs for 
<7i, g 2 , and g x . Including the one-loop exact results, the 
particular form of the gt interaction then implies that the 
RGEs (to all orders) have the following form, 



1 



dgi 
dl 

dm 

dl 

dg x _ 

dl ir(vx + v 2 ) 



^ = -^[1 + %/')] 



2-KVl 

1 



g 2 t [l + 0{g a /t)} 



dgt 
dl 



9t 

IT 



vi + v 2 lv\ 2v 2 



(10) 



The energy (temperature) scale is related to I by T ~ te~ l 
and 0(g™/t n ) denotes higher order terms in the cou- 
plings, g™/t n (a = 1,2, x,t). The plus and minus signs 
result from particle-hole respectively particle-particle di- 
agrams. The one-loop RGEs have been derived in Refs. 
[g4],[14]], but the authors have not noted the particular 
form to all orders. 

The exact solution of (|ITJ) to one-loop order is given in 
Appendix A. We obtain that for comparable ratios of the 
velocities, 1/7 < i'i/v 2 < 7, all couplings diverge at the 
same scale, while for other ratios, gt scales to zero and 
gi, g 2 and g x stay of the order of U. The stability of the 
gt = fixed point follows from the particular form of the 
RGEs to all orders (|l0|). Since all couplings (to all or- 
ders) are multiplied at least once with g t , this fixed point 
is stable for g t — > (any higher order term of the form 
9i 2 x would drive the system to another fixed point). 



IV. BOSONIZATION AND PHASES 

Using bosonization and the above RG results, we de- 
rive the low-energy phases (for an overview, see Fig. 1). 
When doping the half filled ladder, a SC phase with inter- 
chain pairing arises (up to v\jv 2 < 7), see section IV. A. 
In section IV. B, we discuss the g t = phase, present 
well away from half filling, v\/v 2 > 7, and exhibiting 
two gapless modes. In section IV. C, we finally study the 
transitions and crossovers between the different phases. 

Bosonizing the Hamiltonian H, we obtain (see Ap- 
pendix B), 






J = l,2 
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2tt 



[0**10**2 - n^] 



9i 



(2-Kaf 



4tt(0i-0 2 ) . (11) 



A flow to strong coupling of gt results (classically) in a 
"pinning" of B\ — 6 2 = in order to minimize the energy, 
and a single gapless mode, while for g t — ► 0, two gapless 
modes are present. 



A. The interchain-pairing SC phase 

We show that as a result of a finite interchain hop- 
ping t±, doping the half filled ladder, we obtain a phase, 
where SC correlations dominate. The pairing takes place 
between left (right) going particles in chain one and right 
(left) going particles in chain two. For t± -C t, we recover 
previous results p^j. 

When doping the half filled ladder (up to a ratio 
fi/f2 < 7), the couplings grow and eventually diverge. 
However, one should not overinterprct this divergence. 
For small U ', the combination 61—82 in ( |ll|) can fluctuate 
at high temperatures, but for growing g t (i.e., decreasing 
temperature), these fluctuations become suppressed at 
the scale where gt(T c ) ~ t. The 8\ — 82 — then leads 
to a coherence between the two bands and is interpreted 
as a crossover to a new phase. Since T c ~ te~ lc <C t is a 
very low temperature (l c ex t/U), the change of the cou- 
plings as a function of temperature for T < T c is weak as 
long as we keep away from the point of divergence, i.e., 
within the validity of the one- loop RGEs. It is then a 
common practice to set gt ~ t as a strong coupling value 
(and similarly, for the other couplings). 

Using the canonical transformation <I>± = ($1 ± 
$ 2 )/\/2 and n± = (IT ± n 2 )/\/2, the Hamiltonian @ 
takes the form 

H = Hb + Hsg + -fornix, (12) 

where Hb is the Hamiltonian of a massless boson, 



(13) 



H B = Jdx^± 


^-(S^ + jr+n* 


and -ffsc is the sine-Gordon Hamiltonian, 


H S g = J dx I — : 






fft 


V^tf- 


1 


(27ra) 2 C ° u 


J 


and finally, JTmix is a ] 


nixing term, 







(14) 



flmix = / rfx («ift,* + fl x *_ + ^.n+n_) . (15) 

The velocities u±, v ^, and the LLP K± are given by 



u± 



Vi +v 2 



gi + g 2 ± gga 

47T 



c.p 



Dl 



_W2 ± 51 



.92 



4tt 



and 



A"4 



' 27r(fi + i> 2 ) - (fli + 92 ± 2ff x ) 
27r(wi -I- v 2 ) + (51 + g 2 ± 25 x ) 



(16) 



(17) 



(18) 



The mixing term hinders for v\ — v% 7^ a (simple) ana- 
lytic solution of the (classical) equations of motion. Since 
the 0- field is pinned, the current density takes the form 
j = u + K + U + . At half filling the <f> + field is also pinned 
resulting in j = and an insulating phase [g4| . 

Next, we discuss the correlation functions. For repul- 
sive interactions, the charge density- and SC pairing fluc- 
tuations with the most divergent susceptibilities are the 
following ones: The CDW (at a wavevector kp\ + kp2, 
for a comparison, see section II) is given by the operator 



O 



CDW 



d R1 d L1 



d R2 d L2 



*k*« 



*l,* 



R1*M 



tx exp ( iV2n^+ ) cos ( \/2tt8- J 



(19) 



and the SC fluctuations by 



OsC = d R1 dL2 + dR2dLl = *fil*il - ^R2^L2 

ex exp [zV27r6»+] cos (V2n&-) , (20) 

where the rf/j/Lj are the annihilation operators for the 
fermions in chain j Q. The operator Ocdw represents 
an antisymmetric charge density wave and Osc super- 
conductivity with interchain pairing — previously called 
d-wave like due to the antisymmetry with respect to the 
bonding- and antibonding band |26|]. However, the op- 
erator Osc bas odd parity, Osc(~ $) = ~Osc(x), which 
one associates rather with p-wave like superconductivity 
in each band (see Appendix C.l). 

We obtain for the charge density correlation function 
(for a derivation, see Appendix C), 



°cdw( 2; ) cdw 



(0) 



(X X 



(21) 



and for the SC pairing correlation function 

'Ol c (x)0 sc (0)) ex z-V-v, (22) 



where the exponent is 



K, 



7 = 



K + K- 

2u+U- 



( V ty 



(23) 



As a result, the finite t± oc v c _ leads to 7 > 1, implying 
that superconducting pairing correlations dominate. 

The increase of 7 can be understood as follows. The 
pinning of 6- allows us to set n_ = 0. The only coupling 
is then between the $± fields. The $_ field fluctuates 
strongly and affects the $+ correlation function with ad- 
ditional fluctuations thus increasing 7 and stabilizing the 
superconductivity. 

We furthermore find that in the low-energy regime, 
T — > (I — > 00), the LLP K + becomes bigger than one. 
In detail, whether K + is bigger or smaller than one de- 
pends on the following sum (see (fjl§)) 



.91 + .92 + 2g x 



2v\ V 2 



U 



{vi ~ V2) 2 

2v\V 2 



(24) 



Since g x increases for decreasing temperature, the above 
sum becomes negative and K+ > 1. 
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FIG. 2. The exponent 7 of the charge density correlation 
function, oc x ', is shown for U/t — 0.2 and different ratios 
of t±/t. While for small ratios t± <C t, CDW dominate for 
almost all physical values of the doping S, comparable ratios 
t± ~ t favor superconductivity already at lower doping. The 
SC correlations are strongest, when approaching the transi- 
tion to the mixed SC+CDW phase, where the chemical poten- 
tial is close to the bottom of the antibonding band, (7 — > 1.5 
for iii/f2 — > 7). 

It is instructive to expand (E3h leading order in g x /t 
and t±/t (we rewrite g\ and 172 i n terms of g x and neglect 
£/<£), 



7 



1 



1 fvi-v 2 



2 V vi + v 2 



1 , Vi + v 2 

1 +o 9x 

ZirviV2 



(25) 



Neglecting the g x ~ t due to the relatively small prefac- 
tor, we rewrite (|25|) in terms of physical quantities, 



7 



1 



1 ft 



cot(7rn) 2 



8 V * / sin(Trn) 2 - (i_i_/2£) 2 



(26) 



Leading order in the doping 8 = 0.5 — n, we then find the 
particular simple form, 



7 



1 



t{¥* : 



(27) 



The S 2 and the t\ reflect the v\ <->■ v 2 symmetry. Com- 
paring (E71) with the usual exponent if of a LL, we deduce 
that the t± term produces an effective attraction between 
particles of chain one and chain two. 

The exponent 7 as a function of the doping 6 is shown 
in Fig. 2 for different ratios t±/t = 0.1, ... ,1.5 (we have 
set g x — t) [B3[. It becomes smaller than one when ap- 
proaching half filling, 7 rs K + < 1 (a precursor of the in- 
sulating state), and 1.5 near the transition to the mixed 



phase. The insulating state at half filling thus evolves 
(over a phase with weak CDW correlations; this depends 
on the choice of g x ~ t ) to a SC phase with strong corre- 
lations, where the chemical potential lies near the bottom 
of the antibonding band (vi/v 2 ~ 7) Q. A comparison 
of Fig. 2 with Fig. 1 shows that the SC correlations are 
largest where the crossover temperature is highest. 

For small tj_/t we recover previous results: the 
bosonized Hamiltonian goes over to the one studied in 
Refs. p2 26 J and the exponent 7 stays smaller than one, 
7 « K + < 1, and charge density correlations dominate 
(for almost all physical values of the doping S) . 



B. The mixed SC+CDW phase 

For V\jv% > 7, the fixed points of the couplings are 
such that gt — 0, and g\, g 2 , and g x are of the order 
of U (se e Appendix A). We diagonalize the Hamilto- 
nian (|llj) for g t — using a current representation of 
the fields (for a comparison, see Appendix B). Defining 
J = (Jri, Jli, Jr2, Jli) and 



M 



( v\ gi/2n g x /2ir ^ 

gi/27r vi g x /2n 

g x /2ir v 2 g 2 /2n 

\ g x /2ir g 2 /2ir v 2 J 



(28) 



the Hamiltonian (|ll]) can be written as H — tt j dxJ T MJ 
and the new LLP and velocities are determined by the 
eigenvalues of M, which are 



vi +V2 91 +92 , -D+ 
2 4tt 2 ' 

vi +V2 _ £1 +ff2 , D__ 
2 Air 2 ' 



where 



D± = \ [vi - v 2 ± 



gi -92 

2tt 



2tt 



(29) 



(30) 



Since v± — v 2 S> g x , the coupling g x does not contribute 
to D± leading order in U. The velocities and the LLP 
have thus the same form as in two decoupled LL, 



(2^ and K * 



' 2ixvi 



Hi 



2lTVj 



Hi 



93 

2ttvj 



(31) 



where in our case the gj are scale dependent. The new 
basis is "almost" the old one, Qi w III + e^, where 
e ~ U/t (and similarly for H 2 and the fields $j) and the 
current density takes the usual form, 



j = uxKiIii + u 2 K 2 U 2 



(32) 



where the IX,- obey the equations of motion <9 t 2 IL, — 
u^dlflj. It should be noted that here the Drude coef- 
ficient UjKj is not equal to that of the free Fermi gas 

In a short range of ratios, 7 < V\jvi < 8, the low tem- 
perature fixed point of the couplings is such that gi < 
and g 2 > (for v\/v 2 — 7, we find g± = —0.57U and 
.92 = 0.78J7), implying K\ > 1 and K 2 < 1 — a coexis- 
tence of pairing correlations in the bonding- with charge 
density correlations in the antibonding band. Rewriting 
the SC pairing operator Osci = ^ri^li in terms of the 
chain operators dji/hj shows, that there are as well intra- 
and znier-chain SC pairs. We interpret this coexistence 
as a precursor of the SC phase at v\/v 2 < 7, i.e., pair- 
ing correlations first become established in the bonding 
band (at a momentum k^ = 0), and then also in the an- 
tibonding (at a momentum k± = 7r) and phase coherence 
between the bands sets in for increasing g t . In Fig. 1, the 
dashed line shows, where K\ > 1. 

For ratios Ui/«2 > 8, we find gj > and in both 
bands charge density fluctuations dominate, Kj < 1 (for 
V1/V2 3> 8, the fixed points are g\ w g 2 w U). This 
corresponds to the usual ID metallic (LL) phase. 



C. Phase Transitions and Crossovers 

Next, we discuss the transitions (at T = 0) and the 
crossovers between the different phases as a function of 
doping and temperature (see Fig. 1). The transitions are 
accompanied by non-analyticities in the compressibility 
k respectively in the coefficient A of the specific heat, 
C = XT. 

We interpret the energy-scale T c ~ te~ lc (l c oc t/U) at 
which the couplings become of the order of t as the tem- 
perature, where the crossover to the interchain SC phase 
takes place, see Fig. 1 (the high temperature cutoff To 
is fixed but arbitrary). The crossover temperature T c 
increases upon doping, reaching a maximum at a finite, 
t±/t dependent doping, and then decreases to zero when 
approaching the transition point to the mixed SC+CDW 
phase, vi/v 2 ~ 7. 

The compressibility k and the coefficient A of the spe- 
cific heat in the mixed phase have the same form as in 
two decoupled LL p8|, e.g., 



A 

aT 


■1' 


\Ul u 2 ) 


— i| 


f 9i ^ 
\2ttvi / 


)*U> 


( 92 y 

(33) 



where Ao is the specific heat coefficient of the noninter- 
acting system. In the interchain SC phase, the main 
dependencies of k and A come from K + and u + , but 
similarly as for the correlation functions, there are cor- 
rections in v c _. In any case, the values of A and n "jump" 
at v\/v 2 = 7. 



The SC+CDW phase then evolves into a CDW+CDW 
phase, where the crossover takes place at v\jv 2 ~ 8 
(at T = 0). The band edge of the antibonding band 
(i>i/i>2 — ► 00) cannot be treated by the present method, 
since the dispersion relation is there quadratic, ex k 2 . 
However, choosing a sufficiently small U, we can ap- 
proach this point arbitrarily close. We thus conjecture 
that the CDW+CDW phase can be extended up to this 
band edge. 

When only one band is occupied, i.e., v 2 — 0, the phase 
consists of a single LL (with CDW excitations); at the 
band edge of the antibonding band, we then have a tran- 
sition from a one-band CDW to a two-band CDW+CDW 
phase, where for U/t — ► 0, the coefficient A and the com- 
pressibility k diverge as A, k ~ 1/&F2- 



V. CONCLUSIONS 

We have shown that a finite interchain hopping t± has 
the effect of driving the weakly interacting spinless two- 
leg ladder to a SC phase, where left (right) going particles 
in chain one are paired with right (left) going particles 
in chain two. The SC correlations are largest when the 
chemical potential is close to the bottom of the antibond- 
ing band. Between the SC phase at lower- and the CDW 
phase at higher doping, we have found a (new) phase 
with coexistence of SC and CDW excitations. Two phase 
transitions take place (for T — * 0): one from the inter- 
chain SC phase to the mixed SC+CDW phase and one 
from the CDW+CDW to the (one band) CDW phase. 

On the one hand, SC correlations are enhanced by the 
velocity difference v± — v 2 = 2tj_/tan(7rn), but one the 
other hand, this difference suppresses the coherence be- 
tween the two bands, when the chemical potential is too 
close to the bottom of the antibonding band. We ar- 
gue that this sort of competition also takes place in the 
spin-1/2 case 0,01 and generally for (possible) super- 
conductivity in iV-band models. 
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APPENDIX A: RENORMALIZATION GROUP 

The renormalization group (RG) method is a con- 
trolled way of subsequently eliminating (integrating out) 
high-energy modes in a given Hamiltonian. While the 
nonintcracting part Hq is (usually) at a RG fixed point, 
couplings of the interacting part -Hint may grow or de- 
crease under a RG transformation (i.e., when lowering 
the energy). 



A subsequent (perturbative) elimination of high- 
energy modes in H\ nt then results in RG equations, which 
give the change of the coupling when lowering the energy, 
see Ref. [§7|. 

Next, we solve analytically the one-loop RGEs given 
in equation (JIOJ). Keeping only the terms quadratic in 
the coupling constants, we transform the set of the four 
differential equations ([LOj) into one differential equation 
for g x , 



1 dg x 
B dl 



where 



B 



C 



{g x - cuy + du 2 , 



4vi«2 + (vi +v 2 ) 2 
2irviv 2 {vi + v 2 ) 

2( Vl +v 2 ) 2 

^V\V 2 + {Vl + V2) 2 ' 



and 



D = 



-vf + Qvfv2 + 6vfv 2 + Qv\v 
[4viv 2 + (vi + v 2 ) 2 } 



(Al) 



(A2) 



(A3) 



(A4) 



The solution of equation (Al) is qualitatively different 
for D < and D > 0. For D > 0, all couplings diverge, 
while for D < 0, g t scales to zero and the others remain 
of the order of U. 
Solving the equation 



Qx 6 



6x 2 - 6x + 1 = 0, 



(A5) 



for x — V\/v2i we obtain the exact transition ratio v\/v2- 
For comparable velocities, 



1/7 sw 0.14327 . . . < vx/v 2 < 6.9798 . . , 
we find D > resulting in 



g x {l) = U {C + Vd 'tan 



BvDUl — arctan 



(7-1 



(A7) 



and 



gt(l) = Uy/ir(vi + V2)BD 



x < 1 + tan 



BvDUl — arctan 



m 



1/2 



(A8) 



For ratios t>i/t>2 > 7(< 1/7), we find gt — > and 

g x ^{C-V=D)U, (A9) 

for I — > 00. Similar as g x , the coupling constants gi and 
52 stay of the order of U (e.g., U < g x < 1.417). Since gt 
flows to zero and the form of the RGEs to all orders (llfl) 



is such that couplings whatever the order is, are always 
multiplied at least once with g t , the fixed point is stable, 
dg a /dl — > (a = 1, 2, x, t) for I — > 00. For finite U < Vj, 
higher order terms lead to a change of the transition ratio 
Vi/v 2 . 

The one-loop RGEs © have been used in Ref. to 
treat spin-1/2 fermions and the author has noted the two 
different flow regimes, but he has not given the solutions 
of the RGEs. 



APPENDIX B: BOSONIZATION 



Bosonization is a method of rewriting Dirac-fermion 
operators $ in terms of bosonic fields <& and II satisfying 
the commutation relation 



mx),U(v)}=^(x-y). 



(Bl) 



For fermions on a chain or ladder, bosonization applies 
in the continuum limit. A "complete" treatment of 
bosonization is given in Refs. |2|,|28| . 

For convenience, we introduce the dual field of $, 



6{x) 



dx'n{x'). 



(B2) 



For spinless fermions, the bosonization scheme is then 
as follows. We first decompose the fermionic operator \I/ 
into right- and left-movers ^>r and ^l, 



*(x) 



AkpX 



^ R (x)+e-^ x ^ L (x) 



(B3) 



7, (A6) In terms of bosonic operators, the &r/l take the form 



%R/l(x) = 



Vr/l 

V27ra 



exp [ty/^iT^ix) + 9(x))] , (B4) 



where a is a cutoff parameter of the order of the lattice 
constant and the t?_r/l are Majorana ("real") fermionic 
operators (usually called "Klein factors" ) , which are nec- 
essary to fulfill the anticommutation relations of the 
^r/l fields. Currents Jr/l — ^ r/l^ 'r/l trien become 

Jl + Jr = -^d x <5> and J L -J R = -=U. (B5) 

The generalization to N species of spinless fermions is 
straightforward. 

Fourier transforming the Hamiltonian H of the spin- 
less two-leg ladder (two-band model) back to the x space 
and using the above rules, we simply obtain the bosonized 
Hamiltonian (|Tl"|). 



APPENDIX C: CORRELATION FUNCTIONS 

For the spinless two-leg ladder in the bosonized repre- 
sentation (|12|), we calculate (equal time) charge density 
fluctuations 



(0 



,i(3<S> + (x) 



( 



,i0* + (x) 



» 



3 -z/3$+(0) _ / p -i0* + (O) 



(e~ 



'») 



( 



,i/3* + (x),,-i/3* + (0) 



>-!<• 



3 i/3$ + (a:)\ 



(CI) 



and similarly, fluctuations of the SC pairing operator pa- 
rameterized by the #+ field (for simplicity, we drop the 
pinned 6- field) . We make use of the path integral for- 
malism, where expectations values of operators are calcu- 
lated by integration over all field configurations weighted 
with the exponential of the classical action. 

We first revisit the correlation functions of the free 
massless boson and of the SG- model (see e.g. pq|). 



1. Massless boson and SG-model 

The nature of the phases strongly depends on the long- 
range behavior of charge- and SC pairing correlation 
functions. In a single chain, such correlation functions 
are determined by a single (interaction dependent) pa- 
rameter, usually denoted by K (Luttinger liquid param- 
eter, LLP). 

The bosonized action of a single chain of spinless 
fermions is in the continuum limit the one of a massless 
boson, i.e., it is gaussian in the field $ allowing for an 
analytic calculation of correlation functions. It is con- 
venient to perform a Wick rotation to imaginary time, 
r = it. The action then reads, 



S =2K 



l — f dxdT§(-d 2 T + vdl\§, 



(C2) 



where K is the interaction dependent LLP and v the 
Fermi velocity. The Green's function G satisfying, 



1 



is given by 



d^ + vdl )G{x,t) = 6(x)8(t), 



R 2 



G(x,r) = -In 



x 2 + v 2 t 2 + a 



(C3) 



(C4) 



where a is a short distance cutoff and R the radius of 
the integration boundary in the complex plane (we fi- 
nally take R — > 00, corresponding to the usual thermody- 
namic limit). The one-point correlation function is then 
equal zero (e^*^)) = 0. Defining O cdw = & R $> L , the 
long-range behavior of the (equal time) charge density 
correlation function at 2kp is 



°CD\v( ;r ) O CDw(0) 



(X X 



.-2K 



(C5) 



Here, 1 = Air (see above). The SC pairing operator is 
Osc = ^r^l and its correlation function 



Oi c (x)0 sc (0) 



OC X 



-2/K 



(C6) 



A K > 1 therefore implies dominant pairing-, while 
a K < 1 results in dominant charge density correla- 
tions. Note that the SC pairing operator has odd parity, 
Osc(%) = ~Osc{—%) (for spin-1/2 fermions, the parity 
is even). For chiral fermions, the parity transformation is 
^l(x) — » ±Wfl(x) for periodic (antiperiodic) boundary 
conditions. 

The action of the SG-model has the form 



S = ~ f dxdr® (d 2 + dl) $ 



£cos(/3 $) 



(C7) 



where we take /3q < 8tt. The "pinning" term gcos([3o§) 
renders the one-point correlation function of the $ field a 
constant, (e I/3 *^^)=const^ 0. The two-point correlation 
function becomes a constant at large distances (cluster 
decomposition principle) , 



e ip*(x) e -ii39(o)\ 



,i/8*(a0 



-i/3*(0) 



const, 
(C8) 



for x — * 00. The one-point correlation functions of the 
(unpinned) dual field 6 is equal zero and the two-point 
correlation function decays exponentially. 



2. The spinless two- leg ladder 

Next, we calculate the correlation functions of the spin- 
less two-leg ladder making use of the above results. In 
our case, the 9- field is pinned. 

The action of the spinless two-leg ladder (resulting 
from the Hamiltonian (|l2|)) is invariant under the shift 
<&+(x) — > $+(x) + c, where c is any constant. Similarly 
as it is the case for the action of the free massless boson 
we then obtain (for a comparison, see Ref. p5||) 



i£ &*+(*,; 



= 0, 



(C9) 



if Y] . (5j 7^ 0. In particular, the one-point correlation 
function of the $+ field is equal zero, (e^^+^^j = 0, 
i.e., the $+ field is indeed a free (unpinned) field. The 
same holds for the dual field + . We conjecture that the 
mixing term (fLq) for V\ 7^ t>2 is an analytic perturbation. 
Using the Green's function (C4), we carry out the in- 
tegration over the fields II + and $+ and obtain for the 
equal time charge density correlation function 



e i/3<M*) e -i/3*+(0)\ x x -2K+^/4ir I £ S* 



<(*) 



SG 



(C10) 



The expression 5*^ depends (nonlocally) on the fields 
<f>_, II vanishing for v± = v 2 (for simplicity, we drop in 
the following the short-distance cutoff a), 



1® — —i 

w mix 



'/ 



dz x [v c _dZ 1 $-+i— r , 



x — — In 



4tt 



3 T1 n_ 

(x — x\) 2 + u\t 2 n 



dzidz 2 uia^i-Cl) 



K+u+ 



a T1 n_(i) 



x^ln[^ 2 +4 Tl 2 2 ](1^2), (Cll) 

where (1) = (xi,Ti), X\i = X\ — x 2 , T\i = n — r 2 , and 
dzi = dx\d,T\. The average is taken with the SG action 
resulting from (|1J). Since the <&_ field is unpinned, in- 
tegration over $_ gives corrections ex (vi) 2 (and higher 
order) to the exponent of x. 

Next, we carry out a canonical transformation from 
($_,I1_) to (l»_,n_) = (9-,d x $-) (the transformation 
is canonical, since it preserves the commutation relations 
and the integration measure is not changed, because the 
Jacobian of the transformation has determinant one). 
Rewriting 5* ix in terms of the new fields then allows 
us to carry out the integration over Il_. The II_ part 
has the form, 



" / , -,-, I ,-, K + K- , „ N 2 



2K. 



d^n_(i)^n_(i) 



4ttu~ 



(vty 



x f dz 2 tl- {2)d Xl d X2 In [x 2 12 + u 2 + t 2 2 ] 1 . (C12) 

The inverse of the above operator on II_ can be expanded 
in a power series in {vt_) 2 ■ The linear part in II_ reads 

i [ dz^-il) l ^J- d xi ln 



(X-Xl) 2 +U 2 + T 2 



XT + MirT 



-a T1 $_(i) 



Attu- 



Jdz 2 d Xl In [x 2 12 + u\t 2 2 ] d X2 d T2 $_(2)\ . (C13) 



Carrying out the integration over n_ we obtain the fol- 
lowing contribution to the correlation function in (i>L) 2 , 



K_ f (3K+v c _ 
2uZ V 4tt 



dz\ \ d Xl In 



{x-xi) 2 + u\t{ 



X\ + U" + T{ 



(C14) 



Similarly, the pairing correlation function takes the 
form 



e tf39+(x) e -il3e + (0)\ a x -2i3 2 /(4ttK + ) /gS^Jx) 



SG 



(C16) 

where again S^ iK depends (nonlocally) on the fields $_, 
II_, also vanishing for v\ — v 2 , 



n9 _ P_ 
mix 2ir 



j dz x [v°_dt^ 



arctan 



Xi 



U+Ti 



— arctan 



(— ) 



fd Zl dz 2 (v c _d 2 Xl <S>-{\) 



I v_ 



K + u 4 



-fiun-d' 



k. 



x—\ii[x 2 12 +uI_t' 2 
on 



*+'12. 



(1 



(C17) 



The logarithm in (Cll) and the arcustangent in (C17) 



give after partial integration a similar leading order con- 
tribution. The differ ence comes from the iK+ present 
in ( Cll ) but not in ( C17 ). Here, there are only terms 
ex (i>l) 2 . The combination d x arctaXL()9 x m() does not 
give logarithmic contributions after integration. A simi- 
lar calculation as above leads to a decay tx x~ lp , where 



7p = 



3 2 



2-kK a 



1 



K+K- 
2u+u- 



(v c -Y 



(C18) 



In both cases, the remaining part in the $_ fields is 
either real (and does therefore pin the field) or it is imag- 
inary but multiplied with derivatives of logarithmes bee- 
ing strongly peaked at n = and x± = x or x± = 0, 
resulting in a effective contribution ex i(<t^(x) — $_(0)). 
In both cases, the remaining part then becomes a con- 
stant (^ 0) for large x. 



Comparing (C15) with (CIS) we see that one can ex- 
press both correlation function in terms of a single expo- 
nent, 



K+ 



1 K + K- , . 

1 — o. .. (V_ 



(C19) 



2u+U- 



such that the charge density correlation function de- 
cays ex x~P 7 ' 27r and the pairing correlation function 

ex x^ 13 A 271 "?). In our case, /3 2 = 2tt. 



The integral is equal 87r(lnx)/w+. Including all orders in 
t>i, the correlation function finally decays ex x~ lc , where 



(3 2 K + 



1 



2tt 



2u+u- 



y~) 



2 ' 



(C15) 



[1] J.M. Luttinger, J. Math. Phys. 4, 1154 (1963). 

[2] F.D.M. Haldane, J. Phys. C14, 2585 (1981). 

[3] D. Jerome and H.J. Schulz, Adv. Phys. 31, 299 (1982). 

[4] P.W. Anderson, Science 235, 1196 (1987). 



[5] 
[6] 
[7] 
[8] 

[9] 

[10] 

[11] 

[12] 
[13] 

[14] 
[15] 

[16] 
[17] 
[18] 

[19] 

[20] 



E. Dagotto, J. Riera, and D.J. Scalapino, Phys. Rev. 

B45, 5744 (1992). 

T.M. Rice, S. Gopalan, and M. Sigrist, Europhys. Lett. 

23, 445 (1993). 

For a review, see E. Dagotto and T.M. Rice, Sci ence 271, 

618 (1996) and E. Dagotto, |cond-mat/990825C . 

R.M. Noack, S.R. White, and D.J. Scalapino, Phys. Rev. 

Lett.73, 882 (1994). 

C.A. Hayward et al, Phys. Rev. Lett. 75, 926 (1995). 

K. Kuroki, T. Kimura, and H. Aoki, Phys. Rev. B54, 

R15641 (1996). 

M. Uehara et al, J. Phys. Soc. Jpn 65, 2764 (1996). 

E.H. Lieb and F.Y. Wu, Phys. Rev. Lett. 20, 1445 (1968). 

P.A. Bares and G. Blatter, Phys. Rev. Lett. 64, 2567 

(1990). 

M. Fabrizio, Phys. Rev. B48, 15838 (1993). 

D.V. Khveshchenko and T.M. Rice, Phys. Rev. B50, 252 

(1994). 

H.J. Schulz, Phys. Rev. B53, R2959 (1996). 

L. Balents and M. Fisher, Phys. Rev. B53, 12133 (1996). 

H. Lin, L. Balents, and M. Fisher, Phys. Rev. B56, 6569 

(1997). 

H. Lin, L. Balenti 

(1998). 

R. Konik et al. 



and M. Fisher, Phys. Rev. B58, 1794 



9810332 



cond-mat/9806334 and :ond-mat/ 



[21] H.J. Schulz, Phys. Rev. B59, R2472 (1999). 

[22] A. A. Nersesyan, A. Luther, and F.V. Kusmartsev, Phys. 

Lett. A176, 363 (1993). 
[23] The metal-insulator transition in a spinless two-band 



[24 
[25 

[26 



[27 
[28 



[29 
[30 
[31 

[32 

[33 

[34 



model in infinite dimensions has been studied in Q. Si 
et al., Phys. Rev. Lett. 72, 2761 (1994) and L. Craco, 
Phys. Rev. B59, 14837 (1999). 

K.A. Muttalib and V.J. Emery, Phys. Rev. Lett. 57, 1370 
(1986). 

A.O. Gogolin, A. A. Nersesyan, and A.M. Tsvelik, in 
Bosonization and Strongly Correlated Systems (Cam- 
bridge University Press, 1998). 

The effect of disorder on the different phases has been 
studied in E. Orignac and T. Giamarchi, Phys. Rev. B53, 
R10453 (1996) and Phys. Rev. B56, 7167 (1997). 
R. Shankar, Rev. Mod. Phys. 66, 129 (1994). 
H.J. Schulz, in Proceedings of Les Houches Summer 
School LXI, ed. by E. Akkermans et al. (Elsevier, Ams- 
terdam, 1995). 

E.H. Lieb, T. Schultz, and D.C. Mattis, Ann. Phys. 
(N.Y.) 16, 407 (1961). 

M. Yamanaka, M.Oshikawa, and I. Affleck, Phys. Rev. 
Lett. 79, 1110 (1997). 

P. Gagliardini, S. Haas, and T.M. Rice, Phys. Rev. B58, 
9603 (1998). 

The difference of Osc to the bosonized form of Os2 in 
Ref. |E2| is due to an error in J22J (already noted in J26][). 
The difference to the fermionic form of Oscd in Ref. |2q] 
is due to a mistake in |26J (private communication). 
When choosing g x 2> t, we obtain 7 — » 00 and Osc ~ 
const (at T = 0). 

A similar finding has been made in numeric calculations 
(for spinful fermions) in R.M. Noack et al., Phys. Rev. 
B56, 7162 (1997). 



10 



